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A high resolution inelastic neutron scattering (INS) study of fully deuterated Mni2-acetate pro- 
vides the most accurate spin Hamiltonian parameters for this prototype single molecule magnet so 
far. The Mni2-clusters deviate from axial symmetry, a non-zero rhombic term in the model Hamil- 
tonian leading to excellent agreement with observed positions and intensities of the INS peaks. The 
following parameter set provides the best agreement with the experimental data: D — —0.0570(1) 
meV, Bl = -2.78(7) • 10" 6 meV, B\ = -3.2(6) ■ 10" 6 meV and \E\= 6.8(15) ■ 10" 4 meV. Crystal 
dislocations are not the likely cause of the symmetry lowering. Rather, this study lends strong 
support to a recently proposed model, which is based on the presence of several molecular isomers 
with distinct spin Hamiltonian parameters. 

PACS numbers: 75.50.Xx, 75.45.+j, 75.30.Gw, 78.70.Nx 




Single Molecule Magnets (SMM) are crystalline com- 
pounds consisting of nominally identical polynuclear 
molecular complexes of transition metal ions. They ex- 
hibit slow magnetisation relaxation phenomena at cryo- 
genic temperatures. Discrete steps in the magnetisation 
are a typical feature directly related to quantum tunnel- 
ing processes of the magnetisation (QTM). Mni2-acetate, 
[Mn 12 Oi2(OAc)i6(H 2 0) 4 ]-2HOAc-4H 2 0, was the hrst re- 
ported SMM and remains the best studied up till now. 
Since the pionneering work of Sessoli et al. |ll yl Mn 42 - 
acetate has been extensively studied by INS |3J, Ll , EPR 
NMR H , _R aman and Infrared spectroscopy || , 
micro-Hall probe techniques, micro- 
H2I IkH , magnetization relaxation 



specific heat [21 H 
Squid techniques 

01 E3 0j E3 and numerous other bulk measurements. 
Identifying and quantifying the interactions leading to 
slow relaxation and QTM has been the focus of several 
studies in very recent years In this Letter 

we report a high resolution inelastic neutron scattering 
(INS) study on fully deuterated Mni 2 -acetate. In con- 
trast to most other studies no external magnetic field is 
involved. This leads to very accurate values of the rele- 
vant interaction parameters responsible for QTM and al- 
lows us to unambiguously discriminate between the var- 
ious microscopic models which have been pro posed t o 
account for the observed macroscopic behavior 
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[Mni 2 Oi 2 (OAc)_i 6 (H 2 0) 4 ]-2HOAc -4H 2 crystallizes 
in space group 14, and the [Mn 12 Oi 2 (OAc)i 6 (H 2 0) 4 ] 
molecule occupies a position with S 4 point symmetry |2(| . 
The 4 H 2 and 2 HO Ac (acetic acid) solvent molecules 
are incorporated between the Mni 2 -acetate complexes, 
with the 2 HOAc disordered on a fourfold position. Ex- 
change coupling between the Mn 3+ and Mn 4+ ions within 
the complexes leads to a 5 = 10 ground state. A 
strong axial anisotropy with the S 4 axis as the easy axis 



splits the ground state into sublevels from Ms = ±10 
to Ms = and thus creates the energy barrier which 
is responsible for the slow magnetisation relaxation phe- 
nomena at cryogenic temperatures. The appropriate spin 
Hamiltonian to account for this zero field splitting in the 
S 4 point group is given by [ilj : 
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The first term of Eq. is the leading term. Numerous 
parameter sets have been proposed, with those based on 
EPR and INS studies, which are very similar, usually 
considered the most reliable [10. Eq. cannot ac- 
count for all the observed QTM phenomena, in partic- 
ular tunneling through Ms = ±10 in zero field below 2 
K. Additional terms to the ones in Eq. are therefore 
needed. A Hamiltonian including a rhombic term of the 
form 



axial 



E[S. 



(pi $1) 



(2) 



requires a deviation from the crystallographically deter- 
mined S 4 molecular symmetry. This can result from the 
disorder in the solvent structure, and this possibility was 
examined in detail in Refs [lTj, |22( . The presence of six 
different geometrical molecular isomers with slightly dif- 
ferent environments and thus different D and E parame- 
ters was postulated. An alternative model, in which the 
fourfold molecular symmetr y is broken by crystal disloca- 
tions was proposed in Ref. [16j. This model corresponds 
to a broad distribution of site geometries with E values 
broadly distributed around \E\ = 0. Very recent EPR 
and QTM studies strongly favor a discrete distribution 
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FIG. 1: INS spectrum measured on IN5 with an incident 
wavelength of A = 5.9 A at 24 K for 15 h, summed over 
all scattering angles. The peaks are assigned to AAfs = ±1 
transitions, see Figure |3Ji. 



but exact par ameter values to test the specific 
predictions of Refs |l7]> |22( are still missing. 

In an earlier INS study of partially deuterated Mni2- 
acetate the data were of sufficient quality to determine 
the three parameters D, B\ and B\ in Eq. (JTJ 3]. The 
important question of an E term remained open. By an 
upgrade, the time-of-flight instrument IN5 at the ILL in 
Grenoble has gained an order of magnitude in speed, and 
data of significantly higher quality than in Ref. |3j are ob- 
tained. They allow a deeper analysis in terms of Eq. @ 
and lead to a clear distinction between the proposed mod- 
els. 6.5 g of a fully deuterated sample of Mn 12 -acetate 
were used in the present study. The material was pre- 
pared as described in Ref. [2(| using deuterated precur- 
sors and solvents. A cylindrical aluminum container of 
14 mm diameter and 55 mm length was used for the INS 
experiments. Data were corrected for the background 
and detector efficiency using standard procedures. Fig- 
ure n shows an overview INS spectrum (A — 5.9 A) for 
a temperature T — 24 K, at which all the Ms levels of 
the zero field split S = 10 ground state have some pop- 
ulation. The well developed and almost regular pattern 
of inelastic peaks on both the energy loss and energy 
gain sides (positive and negative energy transfer, respec- 
tively) are assigned to AM$ = ±1 transitions between 
adjacent ground state levels. This assignment is straight- 
forward [3j, since only AMs = ±1 transitions are allowed 
in an axially zero field split state, see Figure [Sk- Their 
relative intensities are calculated using Ref. 23] . The 
higher-order B\ and E terms in Eqs. and J5J will mix 
Ms functions, and this becomes relevant in the following. 
The inner part of the INS pattern was measured at the 
same temperature with increased experimental resolution 
(A = 8 A, fwhm of the elastic line 23 /xeV). The result, 
displayed in Figure shows a sharpening of both the 
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FIG. 2: a) INS spectrum measured with an incident wave- 
length of A = 8 A at 24 K for 12.5 h, summed over all scat- 
tering angles. In b) to d) the same data are plotted without 
error bars including theoretical simulations. The same back- 
ground, shown as dash-dotted lines in b), was used in all the 
calculations, b) Eq. (0, D, Bl and B\ values in set I of Table 
c) Eq. parameter set I in Table|I] d) Eq. parameter 
set III in Table □ 



elastic and inelastic features and, in particular, some well 
resolved structure below 0.3 meV on both the gain and 
the loss side. In this spectral range the energy intervals 
and the relative intensities significantly deviate from the 
regular pattern observed in Figure These deviations 
reflect the higher order terms in the spin Hamiltonian, 
and they will now be analysed. The data in Figure 
are reproduced several times without error bars to allow 
comparisons with different theoretical models. The curve 
in Figure [2)3 corresponds to a calculation using Eq. (JTJ 
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TABLE I: Parameter values obtained from fits to the experimental INS peak positions. I and II correspond to the minima in 
the goodness of fit p lot in Figure [I] Set III corresponds to the model proposed in Ref. [r?|. but with D and \E\ parameter 
values from Ref. 22!] . The numbering of isomers is the same as in Refs |17l l22j. The and B\ parameters are the same as in 
set I. 




FIG. 3: Energy level splitting of the S = 10 ground state 
and allowed (energy loss) INS transitions, a) Total splitting, 
purely axial, b) Upper part of a) . c) Splitting of Eq. Q and 
d) splitting of Eq. ||5J with parameters in set I of Table Q 



FIG. 4: Contour plot of the goodness of the fit of Eq. to 
the experimental peak positions as a function of B\ and E 
with constant values for D and B®. The two least-squares 
minima correspond to the parameter sets I and II in Table [I] 



with the D, B® and B\ parameter values in set I of Ta- 
ble [I] This is the best parameter set within the axial 
approximation, and it is the same within experimental 
error to the one derived previously using INS It is 
evident that for energy transfers below 0.3 meV there is 
poor agreement with the experimental data of the present 
study. The upper part of the energy splitting pattern for 
this calculation is shown in Figure|3;. The B\ term splits 
and mixes the Ms = ±2 states in first order. Fitting our 
observed peak positions with the eigenvalues of Eq. J5J , 
i.e. including an E term, leads to the calculated curve in 
Figure 13;. The best parameters are listed in set I of Ta- 
ble [I] The calculated energy pattern is shown in Figure 
The E term splits and mixes the Ms — ±1 states 
in first order, and this splitting is significant. This cal- 
culation shows a remarkably good agreement with the 
experiment, both in terms of peak positions and inten- 
sities. An almost equally good agreement is obtained 
with the parameter set II of Table [I] The two solutions 
I and II correspond to negative and positive values of 



the parameter Bf, respectively. In a purely axial model, 
eigenvalues for equal positive and negative B\ terms are 
degenerate. Figure 0] shows a contour plot, in which the 
goodness of the least-squares fit is plotted in the param- 
eter space spanned by B\ and \E\. The two minima cor- 
responding to the solutions I and II in Tabled are very 
well defined. We confidently conclude that an E term is 
necessary to account for the observed INS data. And in 
contrast to earlier studies we can quantify the \E\ param- 
eter: \E\= 6.8(15) • 1(T 4 meV and \E\= 5.3(9) • 10'* meV 
for solutions I and II, respectively. Of the two parameter 
sets we give preference to set I, i.e. a negative B\ pa- 
rameter, because our value B\ = -3.2(6) • 1(T 6 meV lies 
close to the most recent and ver y ac curate determination 
by EPR of \B\ |= 3.8(1) • 1CT 6 

Of the two models proposed in the literature to ac- 
count for a deviation from axial molecular symmetry we 
definitely rule out the dislocation model |lfj, because it 
corresponds to a broad distribution of sites with E terms 
centered around \E\— 0. This is the energy level pat- 
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tern in Figure |2t and would lead to a broadened version 
of the calculated spectrum in Figure^D. The alternative 
model with a set of six discrete sites for the Mn 12 -acetate 
molecule 

El 

has recently received qualitative support 
from both EPR and tunneling studies HHIi. In Ref. 
[l7j | the D and E parameters of the six isomers were esti- 
mated on the basis of the angular overlap model. Based 
on the same isomer distribution model D and E parame- 
ters obtained by a density functional theory (DFT) calcu- 
lation were recently reported jl^] • The distribution of D 
and E values is very similar to Ref. ^3 1 but the absolute 
values are more reliable. We calculated the energy split- 
ting patterns for the six isomers and the corresponding 
INS spectrum at 24 K, using the D and \E\ values from 
Ref. [22| and our best and B\ values, see set III in 
TableQ] The result is shown in Figure|3i. The agreement 
with the experimental data is remarkably good, consider- 
ing that no adjustment of the D and \E\ parameters was 
attempted. It is not a priori clear why the distribution 
of six species with populations and parameter values as 
given m set III of Tabled gives a discrete INS spectrum 
with well defined peaks. Inspection of the parameter set 
III in Table d reveals that the three isomers 1, 2 cis and 
3 make up 75% of all complexes in the crystal. The iso- 
mers 1 and 3 have very similar D and, more importantly 
in this context, practically identical \E\ parameters to 
our best parameter set I. This, together with some near 
coincidences of transition energies for the other isomers, 
leads to a discrete spectrum which is in good agreement 
with experiment. Our study thus provides support for 
the model proposed in Ref. |T3| and refined in Ref. [22j |. 

Our spectroscopic study thus clearly re- 
veals that the Mni2-acetate clusters in 
[Mni 2 Oi2(OAc)i 6 (H20)4]-2HOAc-4H 2 deviate from 
axial symmetry. A rhombic term in the spin Hamilto- 
nian is essential, and all three parameter sets in Table 
|U reproduce the INS spectra very well. Earlier studies 
based on Landau-Zener (LZ) magnetisation relaxation 
(lif and EPR 0] measurements provided upper limits 
of 8.7 • 1CT 4 mcV and 12.4 • 1(T 4 meV, respectively, for 
the \E\ values in the most strongly rhombically distorted 
complexes in Mni2-acetate. Our study significantly nar- 
rows down the range of possible \E\ values. In contrast 
to LZ and EPR measurements we are not primarily 
probing the complexes with the fastest relaxation or 
the strongest rhombic distortion but the total of all 
complexes. And the \E\= 6.8(15) ■ 10~ 4 meV value of 
set I in Table [I] represents this. If we adopt the isomer 
distribution model, this value accounts for 50% of all the 
complexes; for 37.5% the \E\ value is at least an order 
of magnitude smaller, and for 12.5% it is 13.8 • 1CP 4 
meV. This latter value is in reasonable agreement with 
upper limit estimates of 8.7 • 1CP 4 meV and 12.4 • 1CP 4 
meV from LZ and EPR, respectively. Interestingly, it 
was recently suggested, based on EPR experiments, that 
in deuterated Mni2-acetate the upper limit of \E\ was 



24.8 • 10 -4 meV, twice as high as for the undeuterated 
material |24|. This value is high and outside the range 
of our parameter values. Our experiment was performed 
on a completely deuterated sample, and we took great 
care to handle it in a water free atmosphere . In a 
sample which is not fully deuterated there is additional 
disorder in the acetic acid and water structure, and this 
could possibly lead to a broader distribution of rhombic 
distortions and thus \E\ values. Rhombic distortions 
of individual complexes, likely caused by disorder in 
the solvent structure, are responsible for some of the 
observations in QTM measurements which are not com- 
patible with a purely axial model. The present study 
provides quantitative information about the size of the 
rhombic parameters. This model is not able to account 
for the observation in Mni2-acetate of QTM between 
levels with odd AM values. Possible mechanisms have 
been proposed in Refs 0, 0, H3 • These effects must 
be small and our experiment provides no information to 
quantify them. 
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